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Abstract 

£SJ We model and derive the solution for the problem of a Mode I semi-infinite crack propagating in a 

i-H discrete triangular lattice with bonds having a contrast in stiffness in the principal lattice directions. The 

corresponding Green's kernel is found and from this wave dispersion dependencies are obtained in explicit 
form. An equation of the Wiener-Hopf type is also derived and solved along the crack face, in order 
S_h to compute the stress intensity factor for the semi-infinite crack. The crack stability is analysed via the 

evaluation of the energy release rate for different contrasts in stiffness of the bonds. 

Keywords: inhomogeneous lattice, semi-infinite crack, Wiener-Hopf technique, energy release rate ratio, stress 
intensity factor. 
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1 Introduction 



Models for cracks propagating in lattices have been extensively studied in Slcpyan (2002) . As noted in 



43 

Marder and Gross (1995)] , discrete models enable one to answer fundamental questions about the influence 
of the micro-structure on the crack motion. 

It is possible to consider types of non-uniformity within a lattice, which can bring new effects in the wave 
dispersion and filtering properties of the structure. These non-uniformities may be caused, for example, by 
thermal pre-stress of a constrained lattice whose ligaments have different coefficients of thermal expansion. 
Linearisation near the pre-stressed state may lead to a model of a lattice with contrasting stiffnesses of bonds. 
^sO Such an elastic lattice, containing a moving crack, is analysed in the present paper. 

Problems involving inhomogeneous lattices containing semi-infinite cracks have been recently analysed. Ex- 
amples include those containing particles of contrasting mass Mishuris, et al. (2007)| , and the dynamic extrac- 
tion of a chain within the lattice Mishuris, et al. (2008)] . 



in 

The propagation of the crack may be caused by feeding waves, generated by some remote source, which lead 
to the breaking of subsequent bonds within the lattice. Feeding waves bring energy to the crack-front bonds 
which cause their disintegration one by one and this produces dissipative waves which carry energy away from 
the front. This process has been investigated in |Slepyan (2001a)| , where Mode III crack propagation within a 
. £h square-cell lattice is considered, under the assumption of uniform straight line crack growth, and steady state 

solutions have been obtained. The lattice is assumed to have bonds of identical stiffness which connect particles 
5-j having a common mass. 

The analysis of Mode I and II crack propagation in a uniform discrete triangular lattice was studied in 



Slcp yan (2001b)] . Compared to the case of the square-cell lattice, the equations of motion relate components 



of the vector field of displacements. It has been shown that both the solution to this problem and the explicit 
form of wave dispersion relations can be obtained. 

From these lattice models, it is possible to determine the regions of crack speeds for which we have steady 
state crack motion, and the stability of these states can also be determined by computing the energy release 
rate for the crack. This has been examined in Fineburg and Marder (1999)], [Marder and Gross (1995)] and 



Marder and Liu (1993)] . Another application for a lattice problem is found in |Slepyan et al. (2010)] . Here, the 



model describing a structured interface along a crack with a harmonic feeding wave localised at the faces, was 



used to predict the position of the crack front. Numerical simulations were presented in Mishuris, et al. (2009a) 
Slepyan et al. (2010)] showing that, for a given range of frequencies of the feeding wave, it was possible to have 
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uniform crack growth or, in the non-linear regime of non-steady propagation, to identify an average crack speed, 
which is consistent with the prediction of the linear model linked to the crack propagating steadily. 

The method of solution of these problems involves formulating the discrete lattice problem in terms of 
the Fourier transforms of functions describing the displacements Mishuris, et al. (2009b)] . A Wiener-Hopf 



functional equation is then derived along the crack faces and factorised to obtain the solution. The kernel in 
the Wiener-Hopf equation has the interpretation of the Fourier transform of the derivatives of Green's kernel. 
Similar features also appear in continuum models of cracks which are solved using singular integral equation 
techniques, [add reference] 

In this article, using the method in Slepyan (2001b)| , our main goal is to solve the problem for a discrete tri- 
angular lattice, which is constrained at infinity, containing a semi-infinite crack with additional inhomogeneities 
being brought by the inclined bonds having a contrast in stiffness to the horizontal bonds (see Figure [I]) . This 
contrast in stiffness may arise as an effect of thermal pre-stress, for instance, heating a lattice with bonds that 
have different coefficients of thermal expansion in the principal directions. The problem for a crack in the 
uniform triangular lattice has been studied in [Slepyan (2002) Chapter 12]. 

The plan of the paper is as follows, section [2] includes the description of the problem for a two dimensional 
lattice with a semi-infinite crack, containing particles of common mass with bonds having contrasting stiffnesses 
in the principal directions. In section [3] we rewrite the problem for the lattice using the continuous Fourier 
transform in the crack line direction, and these equations are solved for the displacements of a particle inside 
the lattice. A Wiener-Hopf equation for particles along the line of the crack is then written in section [3] and the 
representation of the kernel function of this equation is derived. An analysis of the dispersion relations obtained 
from the roots and poles of the kernel function is carried out in section [5] In section [6] we solve the Wiener- 
Hopf equation of section [4] In section [7] we evaluate the energy release rate for the crack propagating through 
the inhomogeneous lattice and investigate its sensitivity to the crack speed and to the stiffness contrast. We 
compute the stress intensity factor for the semi-infinite crack propagating through the inhomogeneous lattice 
and show its behaviour, for different stiffness contrasts, as a function of the crack speed in section [8] In section 
[9] we give conclusions on the results presented here. 

In addition to the main text of the article, we also provide Appendices which include more details of the 
derivations of the results contained in main body of the article. In Appendix I, the coefficient of the leading 
order term of the asymptotics for the kernel function near the zero wavenumber is studied. Some comments 
relating to the analysis of the poles for the kernel function of the Wiener-Hopf equation are presented in the 
Appendix II. Properties of the kernel function which are necessary for its factorisation as a product of two 
functions, analytic in different halves of the complex plane, are proved in Appendix III. In Appendix IV, the 
derivation of the solution to the Wiener-Hopf equation is given. Finally, asymptotes in the vicinity of the zero 
wavenumber for the functions obtained after the factorisation of the kernel function are derived in Appendix V. 




Figure 1: The inhomogeneous lattice containing a semi- infinite crack showing the lattice coordinates. The crack 
is located between the rows n = — 1 and n = 0. 
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2 Elastic lattice and the governing equations 



2.1 The notion of Slepyan's crack 

The crack in the infinite triangular lattice considered here (see Figure [I]), is a semi- infinite fault which appears as 
the result of subsequent breakage of bonds in a straight line, and propagates with a constant speed V according 
to j] — x — Vt, (e.g. between rows n = and —1 in Figure [I). Here, x is the horizontal coordinate associated 
with the crack tip node, directed in the line of the propagation of the crack, and t is time. This approach to 
describing the crack propagation, was first used in Slepyan (2001b) . It allows for the reduction of the full 



transient problem to a simple difference problem in the vertical direction following the Fourier transform with 
respect to r), where the forces in the bonds of the lattice can be linked to the displacements and stresses through 
a Wiener-Hopf equation along the crack. Dispersion properties of this lattice configuration can also be obtained 
through the roots and poles of the kernel function involved this equation. The approach is elegant and leads to 
a closed form solution. One of the drawbacks is that the model does not serve the case of small crack velocities, 
as discussed further in the text below. However, for a finite range of values of a subsonic crack speed, one can 
identify a regime of steady crack propagation. In the real physical situation, crack propagation can be treated 
in the averaged sense, and the average crack velocity can be evaluated, as discussed in Mishuris, et al. (2009a)| . 



2.2 Non-uniform lattice 

We consider a triangular lattice containing a semi- infinite crack (see Figure [I]) , where horizontal bonds are 
assumed to have stiffness [i\ and diagonal bonds stiffness //2- In the following, the ratio of bar stiffnesses is 
a = fii/n2- The length of the bonds connecting the particles is normalised to 1. All particles are assumed to 
have identical mass m. 

To characterise the positions of masses within the lattice, it is convenient to use the basis vectors eo = (1, 0) T 
and ei = (1/2, v3/2) T . The in-plane displacement vector at the node x = keo + nei, where k, n are integers, is 
u(i,x). Projected displacements onto the horizontal and vertical axes, associated with the crack tip node, are 
denoted by u and v, respectively. 

We define the constant 
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° R = 2^V2~V 2Ml+Ai2 ~ y^l + (Vi-V2) 2 (2.1) 
which will be determined, in section [5] as the Rayleigh wave speed for a wave that propagates along the crack 
faces aligned with the x-axis, which are traction free. The crack is assumed to propagate steadily through lattice 
with speed V < Cr, and this propagation is caused by the breaking of the bonds between rows n — — 1 and 
n = 0. 

2.3 Normalised crack speed 

Let V* be the dimensionless crack speed defined by V* — V \Jmj ft 2 - Since we consider sub-Rayleigh wave speed, 



using (2.1| we obtain 

y* < g{a) = C±n = 2a + l- v / 3a2 + ( Q -l) ! 

Proposition 1 The function g(a) is increasing for a > 0, with < g(a) < 3/8 and —^(a) < 3/4. 

da 



Proof. For g(a), in (2.2 1, we have 



dg . 1 4a — 1 



da y 2 4^/3a 2 + (a - l) 2 ' 

and with < a < 1/4, it is clear that g(a) is an increasing function on this interval. Now, let a > 1/4. Using 
the inequality 

3a 2 + (a - If > , (2.3) 

leads to 

^(a)>0, for a > 1/4. 
da 
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Therefore, 17(a) is an increasing function for a > 0, and 



Next, owing to (2.3) 



3 

0(0) = , and g(a) < 5(00) = - , a > . 



d 2 9 ( , 1 ( (l-4a) 2 



da 2 4 v /3a 2 + (a - l) 2 V3" 2 + (a - 1) 

which shows -^-(a) is a decreasing function, that takes its maximum as a — > 0: 
aa 

da 4 



□ 



2.4 Dynamic equations for the lattice above the crack 

The equations of motion for a particle inside the lattice (n > 0) are: 

mii = Ml e • {u| x+eo - u| x eo - 2u| x }e + m £ Vj • {u| x+v . + u| x _ v . - 2u| x } Vj . (2.4) 

where vi = ei and V2 = ej — eo. 

Making a change of variable to the moving coordinate system, we set 77 = x — Ft. Then (2.4) in terms of 
projected displacements becomes 

uiV 2 d 2 u 1 

3-9 n ) - 7 [«(»» + V2, n + 1) + «(»? - 1/2, n + 1) + u(r? + 1/2, n - 1) + - 1/2, n - 1)] 

/3 

-•^-[«('7 + 1/2, n + 1) - u(?7- 1/2, n + 1) + v(r) - 1/2, n - 1) - ^(77 + 1/2, n- 1)] 
— a[u(rj + 1, n) + 74(77 — 1, n)] + (1 + 2a)u(?y, 71) = , 

and 

mV 2 d 2 v \/3 

— J2^ n ) ~ + 1/2 ' n + X) ~ " (7? ~ 1/2, 71 + !) - "(^ + V2, n - 1) + u(r) - 1/2, n - 1)] 

3 

~^Hv+ 1/2, n+ 1) +tj(t7- l/2,n+ 1) +w(?7 + l/2,n - 1) + wfa - 1/2, n. - 1)] + 3v(rj, n) = , 

where 71(77,71) and 11(77, ri) are l ne components of the displacement vector u, along row 71 after the change of 
variable. 



3 Solution for the lattice half-plane 



Here, for the intact part of the lattice (71 > 0) we will obtain the solution to the problem of section 2.4 We 
follow the method discussed in | Slcpyan (20 01b)] , which makes use of the continuous Fourier transforms u b and 
v F with respect to the moving coordinate associated with the crack. 

3.1 The continuous Fourier transform of the system 

We take the Fourier transform with respect to 77, and obtain 

[u F {i, 71 + 1) + ti f (£, 71 - 1)] cos(<£/2) - 2[1 + 2a(l - cos(O) + m(e + Wf/^]u F {t n) 

-iV3sm(k/2)[v F (£,n + l)-v F {Z,n~l)] =0, (3.1) 

3[tj f (£, 71 + 1) + tj f (£, n - 1)] cos(£/2) - 2[3 + m(e + i£V) 2 /ti 2 }v F ((, n) 

-i\/3sin(£/2)[M F (£,7i + l)- u F (£,7i- 1)] = 0, (3.2) 
where e is the regularization parameter, < e -C 1. 
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3.2 General solution for the half-plane lattice 

The functions u F and v F are sought in the form 

u F = CA n and v F = DA n 



Insertion of these into (3.1) and (3.2) yield 



C{(A + 1/A) cos(£/2) - 2[1 + 2a(l - cos(O) + m(e + i^) 2 //^]} - h/3sin(£/2)(A - 1/A)D = , 

£>{3(A + 1/A) cos(C/2) - 2(3 + Y)} - h/3sin(£/2)(A - 1/A)C = . 
For non-trivial solutions of C and D, the biquadratic equation in terms of A must be satisfied: 

(A + 1/A) 2 - 2{2 + 2a(l - cos(O) + ^m(e + i£V) 2 }(A + 1/A) cos(£/2) 

3/Lt 2 

+2{1 + cos(0 + 4a(l - cos(O)} + ^m(s + i^) 2 {4 + 2a(l - cos(£)) + — (e + i£V) 2 } = . 

3^2 A*2 

We require that |A| < 1, and this leads to the solutions 



A j =z j ±Jzj-l, j = l,2, 



(3.3) 



(3.4) 



where the sign in front of the square root in A^ has to be chosen so that the condition |Aj| < 1, j = 1,2, for 
(el and V > 0, is satisfied. For j = 1,2 



Zj = [l + 2asin 2 (C/2) + 2 ™(g + W \ cos( ^ /2) 



m 2 (e + i£V) 4 



9 M | 



x asin 2 (£/2) + 



m(e + i£V) 

3/12 



3/i 2 

4sin 2 (£/2) ( nsm 2 {£_/2) 

1/2 



(e + i£Vy 

3/12 



1 - a 



(3.5) 



Then for n > 0, with such choices of A, the Fourier transforms of the displacements have the form 
u F (t,n)\ l / / l ,(A 1 )/ w (A 2 ) /,(A 1 )/,(A 2 ) \ /CiAf 

F (£,n)/ /«( A i)/-( A 2) \-i/ u (A 1 )/„(A 2 ) -iUA 2 )f v (Ai) ) \C 2 A% 

where 



(3.6) 



/„(A) = -\/3sin^/2)(A 2 - 1) and f v (A) = 3(A 2 + 1) cos(£/2) - 2A(3 + m(e + i£F)7/x 2 ) 

3.3 Mode I Symmetry conditions 

For Mode I, using the symmetry relations 

u(r], —n — 1) = u(?7, n) , u(?7, — n — 1) = ~v(rj, n) , 



(3.7) 



we can extend the solution in the upper half-plane (3.6) to the lower half-plane. Then, for n < —1, the Fourier 



transforms of the horizontal and vertical displacements are 

u F (^n)\ 1 / /,(A 1 )/„(A 2 ) / U (A 1 )/„(A 2 )\ /dAr"" 1 

.^(C.njy /«(Ax)/t,(A 2 ) ^i/„( Al )/ w (A 2 ) ifu(^)fv(Ax) ) VC 2 A 2 -«- 

The forces in the bonds directed along — ei and eg — ei, respectively, for a node in the upper half-plane of the 
lattice, have the form 



T(r),n) 



1 v^3 

-[u(rj,n) - u(r)- 1/2, n - 1)] + — [v(r],n) - v(rj - 1/2, n - 1)] 
— [u{rj + 1/2, n — 1) — u(tj, rij\ + -—[v{i], ri) — v(rj -\- 1/2, n — 1)] 



V3 r 



The conditions (3.7) also give 

1 /3 
0) = - [ufa, 0) - uft - 1/2, 0)] + ^— [v(r), 0) + v(r, - 1/2, 0)] , T(r?, 0) - Sfo + 1/2, 0) 



or 



S F (£, 0) = \{l - 0) + ^(1 + e*/ 2 K(£, 0) , T F (£, 0) = 0)e- i «/ 2 



The representation (3.6), then allows S F in (3.10) to be rewritten as 
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S*&0) = 



R J TB(Ai,A 2 )C 



with 



and 



1 

V3 



Ci 
C 2 

B(A l7 A 2 ) = 



R 



(3.8) 
(3.9) 



(3.10) 



(3.11) 



2/„(Ai)/ fl (A 2 ) 

1 - e^ 2 
i(l + e*/ 2 ) 
/„(Ai)/„(A 2 ) / (Ai)/ o (A a ) 
-/«(Ai)/„(A 2 ) -/ M (A 2 )/,(A 1 ) 

4 The Wiener-Hopf equation along the crack 

In this section, we will use the equations of motion for a particle on the upper crack face (n = 0), to derive a 
Wiener-Hopf equation. The kernel function of this equation is also obtained and its properties are investigated. 
This function is used in (a) section [5] for analysing the dispersion properties of the lattice and (b) for computing 
the energy release rate of the propagating crack in section [7j 

4.1 Dynamic equations along the crack face (n = 0) 

Now we consider the equations of motion for particles along n = 0. The equations are similar to those presented 
in section 2.4 However, since particles along the crack face have no bonds corresponding to the vectors —e% 
and e — ei, and are acted on by an external force <f> (see Figure [2]), this should be taken into account. 




<fa) H V + 1/2) 

Figure 2: The particles along the crack face at n = 0, acted on by external forces 4>{rj) and (f>(rj + 1/2). 

G 



Ahead of the crack, for particles along n = 0, the elongation of the bonds corresponding to eo, ei, ei — eo 
and — eo are the same as those for a particle in the intact upper lattice half plane. Also the internal force of 
the bond in the direction of — ei, can be written as 



M(r),Q) =S(r))H(r)), 



M F =S + , 



(4.1) 



where H is the Heaviside function. Using (3.101 and M. , a similar expression for the internal force in the 



bond directed along eo — ei is obtained, M F = S+e lk / 2 . The external forces acting on the particles on the line 



n. = 0, from (3.101, have the form <f>{rj) and <fi(r] + 1/2), which are directed along — ei and eo — ei, respectively. 
Then, the Fourier transform of the equations of motion for the particles along the line n = are 

™( e + i£lOV(£,0) = icos(e/2) U F (^,l)-{2 a (l-cos(C/2)) + l/2K(^,0) 

^2 2 

sm(e/2)^(£, 1) - 1(1 - e-^)(S + + <f> F ) , 



^(e + m 2 v F (k, 0) = sin«/2)«'(e, 1) + \ cos(£/2)z/(£, 1) - \v F & 0) - ^(1 + e-^)(S + + jf) 



V3, 



In order that the expressions (3.6) fulfill these equations, the system 

T F (£, 0) - S F (ti, 0) - -(1 - e-^ 2 )(S + + 4> F ) , T F (£, 0) + S F {t 0) - (1 + e~^ 2 ){S + + ^ 



should be satisfied (see pTj ), (3.2) and (3.10)). Using the representation (3.6) for the displacements when n > 0, 
in the above, leads to 

H(A 1 ,A 2 )C = -A 1 A 2 /„(A 1 )/„(A 2 )R(5+ + <P F ) . 
where the matrix E(Ai, A 2 ) = [H,-j(Ai, A 2 )]f - =1 has the entries 

S 11 (A 1 ,A 2 ) = A 1 / u (A 2 )(/,(A 1 )(cos(C/2)-A 1 ) + v / 3sin(e/2)/„(A 1 )) , 

H 21 (A 1) A 2 )=A 2 A(A 1 )(-N/3/ u (A 1 )(cos^/2)-A 1 )+/ l ,(A 1 )sin(^/2)) , 
5 12 (Ai,A 2 ) =Hu(A 2 ,Ai) and S 22 (Ai, A 2 ) = 3 2 i(A 2 , Ai) . 



Therefore 



C = -A 1 A 2 / t ,(A 1 )/„(A 2 )[H(A 1 , A,)]" 1 ^^ 



Combining this with (3.11) yields 

S F &0) 

with 
and 



AiA 2 



-R T n(A 1 ,A 2 )R(S+ + , 



2det(S(A 1 ,A 2 )) 
n(A X) A 2 ) = TB(Ai, A 2 )adj(H(A 1 , A 2 )) 



adj(H(A 1 ,A 2 )) = 



S 2 i(A 2) Ai) -H u (A 2! Ai)' 



Since 



^ 2 i(A!,A 2 ) S U (A 1) A 2 ) 
det(H(Ai,A 2 )) = Sn(Ai,A 2 )S 2 i(A 2 , AO -Sn(A 2 , Ai)S 2 i(Ai,A 2 ) . 

adj(H(A 2 , Ax)) = -(° 1 J) adj(H(A l7 A 2 )) , B(A 2 ,A X ) = B(A 1( A a ) ( J J 



(4.2) 



(4.3) 



the matrix Jl matrix satisfies the relation 



n(A 1 ,A 2 ) = -0(A 2 ,Ai) 



(4.4) 
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4.2 The Wiener-Hopf equation 

Next we set 



where <S+, S- are functions analytic in the upper a nd lo wer half-planes of the com plex plane £, and as in (4.1 1, 



■M (£, 0) = <S+, whereas the quantity S F satisfies (3.11 ). After inserting this into (4.2), we obtain an equation 
of the form 

L(O5++5_ = (1-L(O)0 F . (4.5) 



Here, 



m = i 



AiA 2 



2det(H(A 1 ,A 2 )) 



R T fi(A 1 ,A 2 )R, 



where Aj, j = 1, 2, are given in (3.4| 



We note that due to (4.4) and the antisymmetry of det(E(Ai, A 2 )) in (4.3) with respect to Ai and A 2 , that 



is symmetric with respect to Ai and A 2 . We put 



A*j = Zj - J zj - 1 ,j = 1,2 , 



and then an equivalent representation for L(£) in terms of z\ and z 2 (see (3.4) and (3.5)) is given by 

3r(A* 1 )r(A*)(z 2 - Zl )y/%=ly/%=l 



H0 = 



where 



r(Al)F(z 2 )^zJ^l - r(A*)F( Zl )^zJ^l ' 
r(A) 



(4.6) 



1, |A| < 1 , 
— 1, otherwise, 



F(z) = 3(cos(£/2) - z f + 6asin 2 (£/2)(l + cos(£/2))(l - z) 
+ — (e + ^Vf[l ~ zcos(£/2) + (1 + cos(C/2))(l - z)] 



The representation for L(£) in (4.6), contains the stiffness contrast parameter a which is found in the function 
F. For the case a = 1, (4.6) is similar to (12.47) of Slepyan (2002) Section 12.3.5], with the only difference 
here being that the functions r(A*), r^A^) have been included to indicate the change of the branch cuts of the 
square roots whenever |Ai2 j > 1. 



In the derivation of (4.6), from p.4| we used 



Aj = Zj -r(A*j)Jz 



together with fact 



1 



+ r(A* j ) y Jz]-l, for j = 1,2 



4.3 Asymptotics of Ai 2 



Here, we obtain estimates for the behaviour of the roots of (3.3 1, in the neighbourhood of zero and infinity, and 
discuss their properties. 

For the auxiliary functions Zj , j — 1 , 2 we have 



3^2 



2cos(£/2) + (-l)V 1 - 4sin^(£/2) 



+ 0(K 1 / 2 ), 3 = 1,2, £^±oc, (4.7) 



where the term in the parenthesis is non-zero for all £. Here, Y q = 0((e + i£^) 2 ) for most values of £, however 
for a periodic sequence of £ = £ m we have this is Y q = 0(1). 
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In the other limit, £ — > 0, we obtain 



where 



^=^^ + 0(|^| 3 ), j = 1,2 



dj = 3(4a - 1) - 16V„ 2 + (-1) J 4M , j = 1, 2 



M = + 3(a - 1) (3a - 4Vf) , 

and V* = Vy/m/ fi2, as above. 

Now we consider the behaviour of Aj 2 at infinity, we have 



A; 



2z, 



0(|^|- 3 ), |e|-yoo 



(4.8) 

(4.9) 
(4.10) 



(4.11) 



that is, Aj — 0(<£ ) as |£| — > 00 and here the sign in front of the square root in (3.4) is chosen to maintain the 
condition |Aj| < 1, j = 1,2, in this limit. 

Near zero _ 

V3£ 



(4.12) 



In order that the correct sign in front o f the above square root is chosen (so that |A 12 | < 1), it is necessary to 
analyse the behaviour of dj, j — 1,2 in (4.9 1 for a > and < V* < Cn(m/ /X2) 1 / 2 . 

The constants c? li2 

Here we show that di 2 do not belong to the negative real axis in the complex plane. Otherwise 



e 

12' 



1 +7^\ d j\ > 1 



and then higher order terms would be needed in (4.12| to determine if | A1.2 1 < 1 as £ — > 0. The proposition 
below is proved for di and the result extends to di. 

Proposition 2 For a > and < V* < C_R(m/^ 2 ) 1 ^ 2 , we Ziaue d\ £ C\K_. 

Proof. We begin by assuming on the contrary that d± is negative, which is equivalent to the inequality 

3(4a - 1) - 16V„ 2 < 4v/4K 4 + 3(a ~ 1) (3a - 4V* 2 ) (4.13) 

with the assumption that the term under the square root in the right-hand side is positive. This gives rise to 
the cases: 



a) 



together with 



or 



b) 



3(4a - 1) - 16V„ 2 > , 
(3(4q - 1) - 16K 2 ) 2 < 16(4K 4 + 3(q - 1) (3a - 4V„ 2 )) 

3(4a - 1) - 16V„ 2 < and 4V, 4 + 3(a - 1) (3a - 4V, 2 ) > 



9 



The Case a). One can obtain the inequality 

The quadratic function of V 2 , in the left-hand side, has the zeros 

3 1 

8 ' 8 



which for a > 0, are both positive. Then (4.13) is only valid if V 2 <E (wf,^!) f° r a ^ 1/4 OT K 2 G { v 2i v D f° r 
a < 1/4. We now show V? lies outside these intervals for a > 0. 
By Proposition [T] 

max{T/ 2 } = 5(a) < u 2 , when a > 1/4 , (4-14) 



where g(a) is given in (2.2 1. 

For a < 1/4, owing again to Proposition [l] 



da 



Thus 



(a)<|, for a>0, g(0)<± and 5 (l/4) = |(l - ^) 
max{V r H 2 } = gt(a) < v\ , when < a < 1/4 . 



V3\ 3 
< - 



This inequality with (4.14| shows that (4.131 is not valid for V* € (0, C r\J m / 1x2} and a > for case a) above. 
The Case b) . The second inequality of this case can be rewritten as 



where now the quadratic function on the left-hand side has zeros 



>0, 



(4.15) 



w 



± - -[a-l±vT=a| • 



When a > 1 the inequality (4.15) is satisfied (since the zeros of the left-hand side are complex), and it 
remains to show that the inequalities 



o Jim , T , 3(4a-l) 
r2 ' -5— and V? > 



v: < 



/'2 



16 



(4.16) 



do not hold. These lead to 



3(4a- 1) 
16 



< g(a) , or \/3a 2 + (a — l) 2 < — — a , 



where the definition of g(a) in (2.2) has been used. From the above inequalities, we have a contradiction to 
(4.16) for a > 7/4. Assume 1 < a < 7/4, then the last inequality leads to 



, 1 11 



This is valid for a e (0, (-1/2 + VS)/2). Hence (4.13) does not hold 



When a < 1, for (4.15) to be satisfied we put 



= - [a — 1 + \fl — a] 



where 



* * — 



, T . 2 3(4a - 1) 
and K 2 > 16 
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Here, as before, the last two inequalities are valid if < a < (—1/2 + v / 3)/2. It remains to check that 



V?>w\ = ha-1 + Vl^] and < — 

2 M2 



(4.17) 



Both of these yield the inequality 



7 

a 2 {a 2 + a - -) > 



Th is is v alid for a > — l/2+ y2 and, since we have assumed that < a < (— l/2+v3)/2, we have a contradiction 
to (4.17). Therefore, (4.13) is not valid for V* € (0, Cn(m/ n 2 ) 1 ^ 2 ] and a > in case 6) above. 



4.3.1 Asymptotes of L(£) 

Here, we obtain the asymptotics of L(£) for £ — » 0. 

Proposition 3 For £ — > 0, we /iawe 

1 



V(£ + i£)(£-i£) 



o(iei) 



-o , 



where 



L 



3 2 



Proof. Using (4.8 1, for £ -> 



2\/3(3 - 8K 2 )(2K 2 ~ 3«)(K 4 - |(2a + 1)K 2 + f ) 



Z2 - Zi 



0(C) , 



and 



F(z j )=B j ? + 0(e), j = l,2 



together with ( |4.6[ ), (4_9J), (|4_10|), ( |4.20| gives ([418}, where 



in 



2y/3(B 2y /d 1 - B iy fc 



3 = 1,2 



Multiplying the numerator and the denominator by -E^v^i + -Biv di , we arrive at (4.19). 
It is shown in Appendix I that L is positive for a > 0, and V < C#. 



□ 



(4.18) 



(4.19) 



(4.20) 



□ 



5 Roots and poles of the kernel function L(£) 

Having obtained the expression for the kernel function M£), we can now derive expressions for the dispersion 
relations of the inhomogeneous lattice. We also show that these relations can be used to predict the behaviour of 
the argument of the function L, which is needed in the computation of the energy release rate for the propagating 
crack in section 7. Let m = 1 and let 



p = 3r(AJ)r(A3)V^-lV«S-l 

and 



so that in (4.6) 



q = r(Al)F { Z2 )^zj^l~ r{kl)F( Zl )^4^l 
z 2 - Z\ 
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We now investigate the zeros of p and q of the above expression. This is carried out by setting e — and 
O = £V, then solving the equations p — and q = for J7. 
The roots of the equation p = are given by: 



/6^|cos(e/4)| 



fif>(0 = [2 M2 si n 2 (^/4) + 4Misin 2 (^/2)] 1 /2 
2^' 



V^6A*2|ain(e/4)| , 

[2 M2 cos 2 (C/4)+4 Ml sin 2 (e/2)] 1 / 2 . 
Next we obtain the roots of the equation q = 0. This equation, after multiplication by 

r(A*)F(z 2 ) + r(Al)F( Zl ) Jz\ 1 , 

leads to 

F{z 2 )\zl-l)-F{ Zl f{zl-l) =Q 

(Z2 - Zl) 



(5.2) 
(5.3) 
(5.4) 
(5.5) 



Solutions of this extended equation are then (|5.4|) , (5.5 1 and 



(D) 



2 (0 



^2^+^-^/ 3^ 2 + ( M i - sin(e/2) | , 
V^|sin(C/2)| , 

2//i + /i 2 + ^ 3^? + (mi - A^) 2 1 sin(£/2) | , 



(5.6) 

(5.7) 
(5.8) 
(5.9) 



and it remains to determine which of the roots of (5.6| are those of (5.1). Note that for each dispersion relation 



Functions (5.4) and (5.5) are roots of equation (5.1) and are in fact removable singularities of L(£). 

The choice of a and £ determines when ( 5.7 H 5.9 1 are solutions of q = 0. Here, we present the results for 



the case a = 1 and give the description of the behaviour of the argument of L. The function fl[ D ^ is a solution 
of ( |5T] ). When 

0.6667 < £/tt < 1.3333 and 2.6667 < £/tt < 3.3333, 



we have f4 ^ s a solution of q = 0. For £/7r satisfying 



0.7614 < £/tt < 1.2386 and 2.7614 < £/tt < 3.2386, 



(5.10) 



(5.11) 



is a solution of q = 0. The inequalities (|5.10|) and (|5.1l|) and their derivations are discussed in the Appendix 



II. 



In Figure |3]^a) we plot the dispersion relations (5.2)-(5.5) and (5.7)-(5.9), as functions of the wavenumber 
£, for j«i = fj,2 = 1. Included in this figure is the ray il — V£ corresponding to the speed V = 0.4504. By 
comparison with the diagram of arg(L(£)) in Figure [3]^b) , we see each intersection of the ray £V with the 
dispersion curves corresponds to a jump in arg(L(£)). 

At approximately £ = tt, we have an additional jump in arg(L(£)). This is a result of the behaviour of the 
factor z 2 — zi contained in L (see (4.6) and (3.4)). This factor, as shown in Figure [3|c) , is either real or purely 
imaginary. As we approach ^ — tt, this factor goes moves from the negative imaginary axis in the complex 
plane, to the positive real axis (note that z 2 — Z\ is the radical function of (3.5) where the positive square root 
is taken). The corresponding effect is a jump in arg(L(£)) of n/2. 

For a lower crack speed, V = 0.21395, we present the dispersion diagram along with the ray t;V in Figure 
|4^a). On the corresponding picture for arg(L(£)) in Figure |4Jb), we see this function jumps when the ray 
intersects the dispersion curves. In particular, the ray t;V intersects the curve corresponding to fi^ in the 
second region of (5.11), which again results in a change in the value of arg(L(£)). Also, in Figure |4]jc) , at 
approximately, £ = 3.47T, we observe that the factor z 2 — z\ moves from the positive real axis to the positive 
imaginary axis. Again, the result we see is a jump of ir/2 in arg(Z(£)) at this point. 

Finally, for \i\ = 1 and /i 2 = 10, (a — 0.1), V = 0.799076 we present the dispersion relations, the plot of 
arg(L(£)) and the imaginary part of the factor z 2 — z 1; in Figure [5|a)-(c). Comparing with Figure[3j we see the 
distance between the dispersion curves along with their height has increased, which can lead to greater absolute 
value of the area under the curve traced by the argument of L. 
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Figure 3: For ^ x = = 1, a = 1, figure (a) shows the dispersion relations ( 5.2 )— ( 5.5 1 and ( 5.7 )—( 5.9 ), along 
with the ray £V for V = 0.4504. The roots of L are given by the solid curves, poles are shown by curves with 



the s mall dashes and the removable singularities are given by the dashed curves. The regions given in (5.10) 

are roots of q 



and (I5.11I) where flj anc ^ ^3^' are r0 °ts of q — are marked by the dashed and solid vertical lines. For 
comparison, we give the plot of (b) arg(L(£)), and (c) the imaginary part of the factor zi — z\ contained in the 



representation (4.6) of L. 
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Figure 4: For ^ = /x 2 = 1, a 
along with the ray £V for V 



= 1, in figure (a) we show the dispersion relations (5.2|-(5.5| and ( 5.7 H 5.9 1, 
0.21395. 



The description of the curves in the dispersion diagram is the same 
as in Figure [3] We again give the plot of (b) arg(L(£)) and (c) the function Im(z2 — Zi), for the purpose of 
comparison. 
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Figure 5: (a) For \i\ = 1, ^2 = 10, a = 0.1, we plot the dispersion relations (5.2 ) — { 5.5 ) and ( 5.7 >— ( 5.9 >, along 
with the ray £V for V = 0.79976. The description of the dispersion curves is given in Figure [3] The plots of 
(b) arg(L(£)) and (c) Im(z 2 — zi) are also presented. 
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6 Solution of the Wiener-Hopf equation 



6.1 Factorisation of the Wiener-Hopf equation (4.5) 



Here, we use the solution presented in Slepyan (2002)| and extend it to the situation when the bonds within 
the lattice have a contrast in stiffness in the principal lattice directions. We define the class J as the set of 
complex valued functions / : R — > C satisfying 



Re(/(0) = Re(/(-0) , M/(0) 
For e > 0, the kernel L £ J and satisfies 

1 



Ind(L(0) 



2tt 



arg(L(+oo)) - arg(L(-oo))] = , 



(6.1) 



arg(L(0) = -arg(L(-0) and L(±oo) = 1 , 



where Ind(L(£)) represents the number of times the contour traced by L(£) in the complex plane winds around 
the origin. The proof of these properties are given in Appendix III. 
Then, the above conditions allow L(£) to be factorised in the form 

L = L<L- . 



where 



with e 



L±(e + i£V,0 = exp(± 
-0, Im £ > for L + and Im £ < for L_ 



~27ri 



lni(e + isV, s) 



ds 



(6.2) 



The function L + (£_) is an analytic function of £ in the upper (lower) half-plane. The introduction of the 
small positive parameter e into L, shifts its zeros and poles (located on the real axis) into the upper or lower 
half of the complex plane. Following this shift, the dispersion curves, discussed in section [5j can be then used 
to identify in which half plane each root and pole is located (see Slepyan Slepyan (2002) | Chapter 2]). 

Suppose £ z and £ p are wavenumbers which correspond to the intersection of fi = l;V with the zeros and 
poles of L, respectively (see (5.2)-(5.5) and (5.7)-(5.9)), on the dispersion diagram. Also at these intersection 



points, let V < V g , where V g — dfl/d^ is the group velocity. Then, in this case £ z and £ p will be located in 
the lower half plane after the introduction of e in L. Similarly, if V > V g then £ z and £ p are shifted into the 
upper half of the complex plane (see |Slepyan (2002) Section 3.3.3]). From the definition (6.2 1, it follows that 
L + (L-) contains all singular and zero points of L which are located in the lower (upper) half of the complex 
plane for e > 0. 



Using (6.1), equation (4.5) can then be written as 



L+S, 



5_ 



(6.3) 



which is the same as (12.58) of [Sle pyan (2002)] . 

Here, the left-hand side of (6.3) is the sum of two analytic functions, one analytic in the upper half of the 



complex plane, the other analytic in the lower half of the complex plane and in order to solve the problem, we 
must separate the terms in the right-hand side in the same way. This is carried out by introducing the load <fi 
in such a way that it allows for this additive split. In Appendix IV, an outline for this procedure is presented 
for the case when the load <j) is chosen so that the right-hand side can be represented as a linear combination of 
Dirac delta functions. 

In what follows, we will consider the situation when the load applied along the crack faces which generates 



term of the form C5(£), in the right-hand side of (6.3) (see Appendix IV for the derivation). Here, C is a constant 
which represents the intensity of the load. Then 5+ and S- have the form 



L. 



C 

me 



5_ = 



e -> +0 , 



(6.4) 
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7 Energy release rate ratio 

Here, we investigate the dependence of energy release rate ratio [Slepyan ( 2002)| on the stiffness contrast pa- 
rameter a. This is defined as 



Go 

~G 



= exp 



arg(i(e + isV, s)) 



ds 



-0 



(7.1) 



where Go is the local energy release rate for the semi-infinite crack propagating through the lattice with speed 
V for sub-Rayleigh wave speed and G is global energy release rate for the crack propagating through the 
corresponding homogenised medium. This ratio describes the dissipation of energy created by breaking bonds 
at the crack front. 

Next we show (7.1) can be obtained from (6.4). However, we note that the derivation of (7.1) does not 
depend on the solution of the Wiener-Hopf equation in the previous section. The definitions of G and Go follow 
Chapter 12] as the energy release rates globally and locally respectively. In particular, as in 
, G can be defined as a limit 



Slepyan 


(2002) 


Slepyan 


(2002) 



G= lim e 2 5_(ifc)5 + (ifc) , 

fe->0 

and the local energy release rate is given as 

G = lim s 2 S-(ik)S+{ik) . 

The global energy release rate G follows from the asymptotes of S± as £ — >• 0. The asymptotes of L± are 



L+(e + i£V,0 
L_(E + ifV,0- 



exp 



== exp 



arg(L(e + isV, s)) 



ds 



for e — > +0, and their derivations are found in Appendix V. 



arg(L(e + isV, s)) 



(7.2) 



ds 



Then, owing to ( |6.4[ ) and ( |7.2[ ). 



c 



where £ — > 0, e — > +0, and then 



</Lo s/l - i| 



= ex P 



arg(i(e + isV, s)) 



ds 



(7.3) 



exp 



arg(L(e + isV, s)) 



ds 



G = C 2 cxp - 



arg(L(e + isV, s)) 



ds 



-0 



(7.4) 



The local energy release rate is then given by the asymptotes of S± when £ — > oo. As a result of (6.1 ) and (6.2 ), 
L± — > 1 as ^ — > oo. In the same limit we have 



5 4 



C 



e-i£ 



C 



^ ->• oo , e ->• +0 , 



and so G =C 2 . Then this and (7.4) imply (7.1). 



7.1 Sensitivity of the energy release rate ratio to the stiffness contrast 



In this section, we set m = fi\ = 1, so that a = 1/^2, and we investigate the behaviour of the energy release rate 



ratio (7.1) as a function of crack speed, when the stiffness ratio is increased (which corresponds to the decrease 
in the stiffness of the inclined bars). The critical crack speed Cr(ol), in this example, has the form (see (2.1)): 



C R (a) 



m—l,fii 



l 

2aV2 



2a + 1 - v/3« 2 + (a - l) 2 



(7.5) 
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Figure 6: (a): The energy release rate ratio as a function of the normalised speed V/Cr(o), based on formula 
( |7.1[ ). (b) The ratio Go/G as a function of the crack speed normalised by the supremum of Cr{o) over a > 
(Cr(0) = \/3/2). We show the behaviour of this ratio for the stiffness contrasts a = 0.05, 0.1, 0.2, 1, and 100. 
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Figure 7: The energy release rate ratio as a function of the normalised speed V/Cr(o), computed using formula 
(7.1|, for the values a = 0.9, 1 and 1.1. 



The lattice Rayleigh speed Cr{o) is decreasing monotonically as a increases. The upper bound of this function 
is Cr(0) = y/3/2, (when the inclined bars are rigid) and for a tending to infinity (the inclined bars are becoming 
much more softer compared to the horizontal bars) the speed Cr(o) tends to zero. In Figure [6]ja), we plot the 
energy release rate ratio as a function of normalised crack speed V/Cr(o), for a = 0.05,0.1,0.2, 1, 100. The 
crack, for low speeds, does not propagate steadily. The intrinsic feature of the model discussed here is that the 
crack propagates uniformly with speed V. If this assumption is not satisfied then the model does not describe 
accurately the propagation of the crack, which is depicted by the non-smooth part of the curves on Figure [6|a), 
e.g. for approximately V/Cr{q) < 0.5 when a = 1. This region of instability is caused by the highly oscillatory 
behaviour of arg(L(£)) for low crack speeds, in (7.1), as observed in section[5j From Figure |6][a), we see the 
percentage of speeds less then Cn(a), where this region is located, is decreasing as a is increasing. 

Contained in Figure (6Tb) , is the plot of the energy release rate ratio against the crack speed which has been 
normalised by Cr(0) = V3/2. For a > 0.05 the ratio Gq/G tends to zero as we approach critical speed Cn(a) 
predicted by (7.5). For a > 0.1, in the region of instability, the ratio Gq/G increases until it reaches a speed 
where there is a global maximum. For speeds greater than this, we observe the ratio follows a smooth curve 
as it decreases to zero at the critical Rayleigh speed (7.5). In this region, this behaviour shows that the model 



describes the steady propagation of the semi-infinite crack through the inhomogeneous lattice. 

In Mishuris, et al. (2009b)| , the dynamic problem for the propagation of a crack situated in a square-cell 
lattice, whose rows contain particles of contrasting mass, is discussed. The sensitivity to the contrast in mass, 
of the ratio of the local energy release rate Go (for the lattice) to the global energy release rate G (for the 
corresponding homogenised material) was also investigated. It was seen the values of Gq/G, for high crack 
speeds, increased monotonically as the contrast in mass was reduced. Here, as shown in Figures [6]j a), (b) and 
[7j we do not have any monotonicity in the behaviour of Gq/G for different a. For example, in Figure |6ja), the 
curve for a — 0.2 intersects those for a — 1, 100 at approximately V/Cr{q) = 0.83 (where we are in the region 
of stability for all these curves). For speeds higher than this value, the ratio of Gq/G for a = 0.2 is greater than 
that for a = 1 and a = 100. A similar example can be seen for the curve corresponding to a = 0.1. 

It was also seen in Mishuris, et al. (2009b) , for low crack speeds corresponding to the instability region of 
the model, the values of Gq/G where not monotonic as a function of the mass contrast parameter. A similar 
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feature is also observed here for low crack speeds as the stiffness contrast parameter a is varied. In Figure 
[7j we plot the energy release rate ratio for a — 0.9, 1 and 1.1. For speeds V/Cn(a) < 0.5 (the region of 
instability), there exists regions where all curves are intersecting one another and overlapping. For instance, 
when 0.35 < V/Cr{oi) < 0.5, again there is no monotonicity in the behaviour of G /G as a varies. We note 
that the critical speed in the scalar case Mishuris, et al. (2009b)] , was a material constant and did not depend 
on the contrast ratio. 



8 Stress intensity factor in the homogenisation approximation 

In this section, we derive the expression for the Mode I stress intensity factor for the semi-infinite crack in the 
case when the load (j> applied along the crack faces generates the term C<5(£) in the right-hand side of (6.3), with 



C being the load intensity. The solution of this problem is given in (6.4) 



Using the fact that the energy required to break the crack front bond Go = C (see the previous section), 



we compute the inverse Fourier transform of iS+ in (7.3), so that to leading order 



1 Gp 

LoTTT] 



exp 



arg(L(e + isV, s)) 



ds 



for rj > 



(8.1) 



which depends on a and here Lq is given in Proposition [3] Here, the expression /i2<S(?7) represents the tensile 
force in the inclined bonds ahead of the crack. Let the normal traction ahead of the crack be a, then 



3Gq 



j ex P , 



arg(i(e + \sV, s)) 



ds 



(8.2) 



The Mode I stress intensity factor Kj, in the far-field is defined by the formula 

Ki = lim \J 2irn a . 



On comparison of (8.2) and (8.3), we have the stress intensity factor K; is 



K I (a) = fi 2 \l -j— exp 



arg(i(e + isV, s)) 



ds 



(8.3) 



(8.4) 



Now, setting fix = 1 and allowing /12 to vary, we investigate the behaviour of Ki as a function of V and 
a. The expression for Kj/^/Gq is given as a function of the normalised speed V/Cn(a) in Figure |8|a) and as 
a function of the normalised speed V/Cr(0) in Figure |8fb) for a = 0.05,0.1,0.2, 1, 100 (where V < G^(a), see 



(7.51) 



As discussed in the previous section, at low speeds the model presented here does not describe the propagation 
of the crack through the inhomogeneous lattice accurately. This is due to the highly oscillatory behaviour of the 
argument of L(£) and it leads to regions of instability, represented by the non-smooth behavi our in Gn /G f or 
low speeds as seen in Figure^ Since the expression Kj/^/Go contains the ratio y / G/Go (see (7.1) and (8.4)), 
we also see a similar behaviour in Ki/y/Go in Figure pMa) and (b). For those speeds in this unstable region 
there is little physical significance of the results presented in the figures below. As we increase the speed of the 
crack, we leave the region of instability, and the stress intensity factor becomes a smooth function of the crack 



speed, and tends to a finite value for V tending to the Rayleigh speed Cn(a) in (7.5), as shown in Figure |8[a). 
For a < 1, in the vicinity of G#(a), Kj is a monotonically decreasing function for increasing V. When a > 1, 
outside the instability region, K] passes through a global minimum before the converging to a finite value at 
Cu(a). From the figure, we can see that overall the behaviour of the stress intensity factor is increasing as we 
decrease the stiffness contrast a. 



9 Conclusions 



The paper has brought together analytical insight on a propagating semi-infinite fault in an elastic lattice with 
the computational experiments and physical interpretation of the fields in such a lattice structure. One of the 
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V/C R (a) 




V/C R (0) 



Figure 8: The normalised stress intensity factor Kj/^/Go as a function of (a) the normalised crack speed 
V/Cn(a) and (b) the normalised crack spee d V/Cr(0), based on formula ( |8.1[ ), for a = 0.05,0.1,0.2,1,100. 
Here Cr(0) = a/3/2 is the supremum of (7.5| for a > 0. 
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important issues discussed here is the crack stability for different regimes of the crack speed and different values 
of elastic parameters of the anisotropic lattice. 

Although it appears that the analytical solution has a serious limitation due to the assumption of the steady 
crack propagation, it has been shown in Mishuris, et al. (2009a)| that in the non-steady regime the averaging 



procedure appears to be viable. In this case, the averaged solution follows the prediction of the linear model 
constructed for the case of the steady crack propagation through the lattice. 

It is also noted that the low-speed unsteady regime, predicted by the model, appears to be consistent with the 
results obtained by other approaches (see, for example, [Marder and Gross (1995)] ), and it has a clear physical 
interpretation of an accelerating crack at the initial stage (low speed) of the crack advance. 

In particular, the crack stability for different values of the lattice contrast is a challenging issue, which has 
been comprehensively addressed here. 
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Appendix I: The sign of Lq 



In this section we show that for V < Cr and a > the constant Lq in (4.191 of Proposition [3j section |4j is 
positive. 
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Note that by Proposition [T] for V 2 < C\m/^2 

3 - 8V? >0 , 2V? - 3a < 



and 



Also 



K 4 -^(2o + i)K 2 + ^>o. 



2a + 1. 



V* + o [a + r > , 



(1-1) 
(1-2) 

(1-3) 



since here the quadratic function of V^ 2 was shown in the proof of Proposition[2]to be positive for V 2 < C\Mj fi 2 . 



Proposition 4 For a > 0, V„ 2 < C\mj [i 2 the constant Lq in (4.19) is positive. 

Proof. Part 1. Consider the cases when a) a > 1, < V? < C\mj [i 2 or b) — 1/2 + \/2 < a < 1 so that 
- 1 + VT^} < K 2 < C 2 R m/v2. 



For the parameter values in both a) and b), Mq defined in (4.101 is nonnegative. 

We show that B\ is negative, which implies B 2 is also negative. The inequality B 1 < gives 

2K 4 - 3(q + 1)K 2 + {j« > V?M • 



(1-4) 



Here, this is valid for 



which is the case, since 



3a + 3 - 3v / a 2 TT 



v: < 



C 2 H m _ 2a + 1 - y/3a 2 + (a - l) 2 3a + 3 - 3%/a 2 + 1 
M2 4 4 

Indeed, the above provides us with 

3Va 2 + 1 - v^a 2 + (a - l) 2 < a + 2 , 
and the left-hand side is positive as a result of 

9(a 2 + l) >3a 2 + (a-l) 2 . 



for a > 



(1-5) 



Thus from (1-5 1 we have 

(a 2 + l)(3a 2 + (a- l) 2 ) - (2a 2 -a + 1) = 2a 3 > , for a>0 



Then (1-4) leads to 



16 



(8^ 2 -3)(4K 2 -3a) 2 <0 



if V? < 3/8, and this holds by Proposition [Tj Therefore Bj < 0, for j = 1, 2. This together with (jT^TJ, (O) and 
( L3| shows L > 0. □ 

Part 2. Consider the cases when a) < a < -1/2 + \/2, < K 2 < C%m/n 2 or 6) -1/2 + y/2 < a < 1, 

< K 2 < §[«- 1 + V 7 !^]- 

For the parameter values in cases a) and b), Mq in (4.10) is purely imaginary and can be written as 

M = iRo , R = y/3(l - a) (3a - 4^ 2 ) - 4 V; 4 , (1-6) 

where according to the proof of Proposition [2j Rq > since V? < C^m/ ' [i 2 . Then 

^2 = <ii and B 2 = B\ , 
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and from (4.19) 



In 



\di\RB(B iy /dD 



V3(3 - 8V 2 )(2V 2 - 3a)(K 4 - §(2a + 1)K 2 + ¥) 



By (I-ll, (1-2) and (1-3) it remains to show that Re(BiVdi) is negative. We have 

Re^v 7 ^) = Re(Si)Re( v / rfi) + Im(Bi)Im(-v/^) , 

with 

Re(-Bi) = -^K 4 + -y-K 2 - , Im(Bi) = > . 



Note that due to (1-6) and (4.9), — ~ < arg(v3i) < 0, therefore 



The roots of Re(i?i) are 



Considering case b), we have 



Re(Vdi) > , Im(Vdi) < . 

- [a + 1 ± \JcP- + 1] > for a > 
V; 2 < u 2 , for 3/4 < a < 1 , 



(1-7) 
(1-8) 



and then Re(£?i) < 0. This with (1-7), (1-8) implies Re(_BiV3i) is negative and part b) is proved. 



Returning to part a), for Re(i?i) < we require 

2a + 1 - ^/3a 2 + (a - l) 2 



< v_ 



which was shown to be valid for a > in the proof of Proposition [2] section [4] Therefore, the proof is complete. 



□ 



Appendix II: Poles of L(£) 



Here, for [ii = f/,2 = 1, we verify the inequalities (5.10) and (5.11 ) of section[5] which indicate when the functions 
(5.8) and (5.9) are solutions of ( |5.1[ ), respectively. 

For f e [0,471-], sign{r(A^)/r(A^)} = 1 when £, 2 V 2 = (Cl^) 2 . By direct substitution of (H into (O) we 



obtain 



where 



1 



l/a , Dh2/e2 -±|isin 2 (e/2)|tan(e/2)||4cos 2 (e/2)- 11^(0-^(0}, 
T&) = sign{2(l + (-l)-'sign{cos(e/2)})cos 2 (e/2) - 1} , j = 1,2 



(II-l) 



and the sign in front of the right-hand side of ( II- 1 ) depends on the sign of r(A*). The factor in front of the 



curly braces in (II-l ), is zero when £ = 0, |tt, |7r, 27r, |tt, ^7r, and 47r, and is singular at £ = it and 37r. 



We now show that (II-l ) is zero in the neighbourhood of these poles, by solving the equation 

Ti(0-T a (0=0. 

For £ = 7r, 37r it is clear this equation is satisfied. When £ ^ tt, 3tt, we are left to determine when the inequality 

1 -4cos 2 (£/2) > , 

is satisfied. This occurs for £ e [ f 7r, f 7r] U [|tt, ^§tt]. The above is confirmed in Figure[9j 
Next, insertion of (5.9) into (5.1) yields 



h/3sin 2 (£/2)|sin(£/2)|(2cos(£/2) - 1 + VS)^V3 - 2 sin 2 (£/2) 



{r(AJ) + r(A*)}. (II-2) 



2W2cos 2 (£/2) 



The factor outside the curly braces in (II-2) has the roots £ = 0,0.7614tt, 1.2386tt, 2tt, 2.7614tt, 3.283tt and 4tt. 



We must have sign{r(A^)/r(A2)} = -1 when V 2 £ 2 = (Q^) 2 - Figure |10| shows the plot of sign{r(A;|;)/r(A|)} 



in this case, and we see that this condition holds provided £ satisfies the inequalities of ( 5.11| ) 
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Appendix III: Asymptotic properties of L(£) 

We now state and prove some properties of the kernel function £(£), which are used in the analysis of the 



Wicncr-Hopf equation (4.5), in section 6 



Lemma 1 L(£) — > 1 as |£| oo. 



Proof. According to (4.11), for |£| — > oo, we have 



A *j = ^-+0(\z j \- 3 ), and r(AJ) = l,j = l,2, 



z?-l = z j + 0(\z j \- L ), (III-l) 



F(^) = 3z 2 - -(e + i^) 2 z,[2cos(^/2) + 1] + 0(|C| 2 ) 



Since 0j = 0(|£| 2 ) for |£| -> oo, then 



□ 



F(z 2 )^/z 2 -l- F( Zl )yJz*-l = 3z lZ2 (z 2 - zi) + 0(|C| 4 ) , 

3(z 2 - z^zl-l^zl-l = 3z lZ2 (z 2 - Zl ) + 0(|C| 2 ) . 
The leading order term in both of the above estimates is 0(|£| 6 ) for |£| —> oo, therefore 

L(0 = i + o(ier 2 ), i^i->oo. 

Next we state 

Lemma 2 For t/ie product of two complex functions /(£), <?(£) G ^7, we have (/<?)(£) G ^7- 
This leads to 

Lemma 3 For £ G M, /or e > 0, L(£) G J. 
Proof. The term 

! 7 l( e 2_ C 2 y 2 +2i ^ ) g j 

since it has an even real part and odd imaginary part with respect to £. This together with Lemma[2]and (3.5) 
shows that Zj G J, j — 1,2. and hence, after a second application of Lemma [2] with (4.6) we complete the 
proof. 

Now we consider the index of the function £(£) which defined by 

Ind(L(0) = i [arg(L(oo)) - arg(L(-oo))] . (III-2) 

This quantity can also be interpreted as the number of times the path traced by L(£) in the complex plane for 
£ G (—00,00) winds around the origin. We have 

Lemma 4 For e > 0, we have Ind(L(£)) = 0. 

Proof. According to section 6, if £ = £* represents a root or a pole of L(£) such that V < V g (V > V g ), 
then the rcgularisation parameter e shifts this root or pole into lower (upper) half of the complex plane. 
Correspondingly, £ = — £* is then a root or a pole of L(£) such that V > V g {V < V g ) and this will be located 
in the upper (lower) half of the complex plane after the regularisation of £(£). Therefore, L(e + i£V,£) has 
an equal number of roots in the upper and lower half planes and an equal number of poles in the upper and 
lower half planes, and these numbers are independent of e. Hence it is sufficient to show that for large e, 
Ind(Z(e + i£V,£)) = 0. 
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First assume £, e ^> 1, so that |e + i£V| 3> 1, then 
Z ; = 0( £ 2 + £V 2 ), A* 

and so 



1 

2z~ 



O 



(e 2 + ev 2 ) 3 



r(A*) = 1 , 3(z 2 - zx) Jz 2 - lJzj - 1 = 3z!Z 2 (z 2 - Zl ) + 0(e z + £ V ) 



and 



F(z 2 )y/4^1-F(z 1 )y/4 



2,z x z 2 {z 2 -zi) (e + i£V) 2 (z 2 - z x ){2 + cos(£/2)) 



0(e 2 



Thus, from (4.6 1, 



m(e + i^) 2 [2 + cos(e/2)] 

3^2^1^2 



( £ 2 + £2y2)2 



2x^2 



(III-3) 



where according to (3.5 1, z\Z 2 is 0((e 2 + £ 2 F 2 ) 2 ). The remainder estimate here is uniform for large e and £. 

The plot of the contour traced by L(£) for large e and f, with m = fX\ = \i 2 = 1 (a = 1), V = 0.3941, 
can be found in Figure 11 In this figure both the plot given by (4.6| and the leading order term of (III-3) are 



presented, and it can be seen there is a good agreement between both computations. When e>^, from (III-3) 

/i 2 (2 + cos(£/2)) 



1 



2i£V 



and this asymptotic expression (III-3) implies that for large e the contour traced by L does not cross the 
negative real axis, since to leading order the above real part is always positive, whereas the imaginary part can 
cross the real axis for £ £ (—00,00). Similar behaviour can be seen in Figure 11 where the contour traced by 
L(e + i£V, £) for £ g (— 400-7T, 4007r) in the complex plane forms a closed loop which passes through 1. Therefore, 
Ind(i(e + i£F, £)) = in the limit e — > 00 and this also holds for all e > 0. □ 



The next two corollaries provide us with additional properties of L(£) needed for the factorisation of (4.5 1. 
As a corollary of Lemma [3] we have 

Corollary 1 For £ E R, e > 0, 

arg(L(0) = -arg(L(-0) . 
Owing to (III-2), Lemma |4] and Corollary [I] we also have: 

Corollary 2 arg(L(£)) = as £ — > ±00. 



Appendix IV: Solution of the Wiener-Hopf equation (6.3) for a par- 
ticular load 



In order to determine the functions S + and 5_, as solutions of the Wiener-Hopf equation (6.3 1, for the sake 



of simplicity, we consider a certain type of "load" (f> F which produces simple examples for the additive split 
P + + P- in the right-hand side of (6.3). 



Let V+ CP-) be the set of poles and Z + (Z-) the set of zeros of Lie + i£V, £) that are located in the upper 



(lower) half plane. The first factor in the right-hand side of (6.3) is singular at the zeros of L_ and the singular 
points of L + . Assume ^ € P_, and 



A r , 



[C_ e -i(£ -£-)]« 



a > 0, e ->• +0 , 



where = l/(Vg — V) for £ — > £ p . Also, suppose for e 



L_(e + ^V,0 ~ A z [C+e - i(£ - £+)] b , 6 > 0, e -> +0 , 



(IV-1) 



(IV-2) 
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0.02 




-0.02 1 1 1 1 1 1 1 1 1 

0.995 1 1.005 1.01 1.015 1.02 1.025 1.03 1.035 

Re(L(e + i ^)) 



Figure 11: The contour traced by L(e + in the complex plane for e = 10, m = fii = /i 2 = 1, (a = 1), 

V = 0.3941, £/ir £ (—400,400). The computations for the solid curve are based on formula (4.6) and those for 
the dashed curve are based on the leading order part of ( III-3 1 . 
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for £ — > £ + and = — V g ). Here, £~ is a pole and is a zero of L located in the lower and upper half 
of the complex plane, respectively, and the constants A p and A z depend on a. 

We consider two examples of the function <fi: 

a) 

4>a(v) = 4>{v) = -C A (2C-£) a exp[(C^£ - i^)v]H(-v) , (IV-3) 
and after taking the Fourier transform, we have 



"A,- 



(2(-e) a C A 







and b) 
so that 



K-e + iK-ST)]' ^ 

>b(v) = Q(V) = C B (2C+e) b exp[-(C + e + i&W (v) 
(2C+e) b C B 



(IV-4) 







[C+e-i(£-6 + )] 

In cases a) and b), we have Ca and Cb are constants representing the intensity of the load. Also 

lim <Pa(v) = nm <Ab(»?) = ■ 



(IV-5) 



(IV-6) 



Right-hand side of (6.3) for the load of case a). We insert (IV-4 1 into (6.3) and taking the limit as e — > +0 
and noting that due to ( IV-6 1 



lim 



= 



we have using (IV-1 ): 



L+S+ + — = AX A lim ( 2 C-g)° 



Here, according to Slepyan (2002) Section 2.2 A], the above limit converges to the Dirac delta function, therefore 



L+S+ + — ~ 2tt A p C A 8{i -£,-) = A P C A 



1 



1 



Right-hand side of (6.3) /or i/ie Zoad o/ case 6). Similarly, using (IV-2), (IV-5) and 



by [Slepyan (2002)1 Section 2.2.4], we have 



lim 



C B 



lim (As +L+ = 

E->-+0 



(2C+e) 



+o A z [Q +£ + i(£ - 6 + )]"[C+e - i(e - if)] 

1 1 



Cb 

A* 



, £ -> . 



General solution of (6.3). Consulting the above cases a) and b), by linear superposition the right-hand side 



of (6.3) has the form 



L+S+ 



2J A p 



for e — >• , where and £? p are arbitrary complex constants. 
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Then the functions S+ and <S_ are given by 



MO 



5_ = L_(0 



g 2 



(IV-7) 



E 



E 



where e +0. In section [6j ( |6.4[ ) follows by considering terms in (IV-7) which correspond p_ 
in p\^2| and (fiV-31. 



and v = 1/2 



Appendix V: Evaluation of L± for £ — > 



Now, using the Cauchy-type integral of (6.2) and the asymptotic representation (4.18) for L when £ — > 0, we 



derive asymptotes of the functions L± near zero, found in ( |7.2[ ). The logarithmic term in (6.2) is rewritten as 

In L{e + i£V, = ln|L(e + %V, 01 +i arg(L(e + i£V, £)) . 

Then 



I±(e + i^V,0 = exp ± 



1 

2tti 



]n\L(e + isV,s)\ 



ds 



xexpl ± 



2tt 



arg(L(e + isV, s)) 



)• 



(V-1) 



The first exponent on the right, by the Cauchy theorem, defines the modulus of L(e + i£V, £), whereas the 
second exponent, owing to the expansion 



1 ^ tn 

/ < „n+l 



s — £ * — ; s 



yields 



exp 



arg(L(e + isV, s)) 



d S )~e X p(±^| 

«/ — ( 



arg(L(e + isV, s)) 



The integrand in the right-hand side is an even function as a result of Corollary [TJ and so 



ds 



exp 



arg(L(e + isV, s)) 



exp 



(4 



arg(L(e + isV, s)) 



f/.S 



Therefore, allowing e — » +0 in (V-1), and using Proposition [3j we obtain (7.2) 
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